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SUMMARY 


An analysis of available airplane measurements of vertical gust 
velocity is presented in order to arrive at a simple description of the 
frequency and intensity of gust velocities experienced by airplanes in 
operations. For the purpose of application? tx®. missile operations, the 
results obtained are modified to eliminate the effects of storm- avoidance 
procedures normally followed in airplane operations. The frequency dis- 
tributions of gust velocity are then converted to a form appropriate for 
use in power spectral response calculations. Methods of applying the ' 
results to the estimation of the large and the small repeated loads in 
missile operations are then developed. Simple methods of estimating the 
gust loadings that will be exceeded with a given probability are pre- 
sented in terms of missile response parameters and turbulence parameters. 
The limitations of the present results are also discussed briefly. 


INTRODUCTION 


The effects of atmospheric turbulence on airplane structural loads 
have been of concern for many years. Recently, it has become increas- 
ingly clear that certain types of missiles and unmanned vehicles are 
also sensitive to turbulence in regard to structural loading and control 
problems. It is the purpose of the present paper to extend recent result 
on the estimation of gust loads for airplane operations (refs. 1 and 2) 
to the case of missile operations. In reference 1, initial descriptions 
of the frequency and intensity of atmospheric turbulence and their vari- 
ation with altitude were derived in terms of discrete or derived gust 
velocities, and methods of applying these data to load calculations for 
airplane operations were presented. More recently the development of 
random-process applications to gust response problems has, in turn, led 
to efforts to utilize these data in order to establish an appropriate 
description of the turbulence environment and a procedure for response 
calculations in terms of the power spectra of turbulence (ref. 2). This 
procedure provides a more realistic representation of the turbulence 
field and furthermore is more suitable to the treatment of missile sta- 
bility and elastic dynamics. 
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In the present paper, use is made of data on atmospheric turbulence 
obtained from airplane operations. The application of airplane gust 
data to the calculation of gust loads on missiles involves a number of 
problems among which the following two are of importance. First, modi- 
fications to the atmospheric- environment data obtained from airplane 
surveys are required in order to account for the effects of the storm- 
avoidance procedures normally followed in airplane operations and not 
applicable to missile operations. Second, and perhaps a more serious 
problem, is that concerning the flight-path angle of the missile. For 
flight paths that are moderately inclined to the horizontal, the indi- 
cations of the approximate isotropy of atmospheric turbulence (refs. 3 
and 4 ) suggest that the airplane data would apply reasonably well. For 
flight paths that are more near vertical, however, serious questions 
exist as to the applicability of gust data obtained from airplanes in 
horizontal flight. However, no adequate alternative appears currently 
available for this vertical-flight case. Thus, the present study might 
be considered to apply best to missile operations in flight paths similar 
to those of airplanes or in -moderately inclined flight paths and to 
apply only in a very crude way to near-vertical flight paths. 

This paper presents the results obtained from an examination of 
available data on the frequency and intensity of atmospheric gust veloc- 
ities and their variation with altitude and, in this respect, brings up 
to date the results reported earlier in reference 1. These data are 
then adjusted for the present purpose of missile application to account 
for the effects of airplane storm-avoidance practices. The distributions 
of gust velocities are then converted into a form appropriate for use in 
power spectral response calculations in accordance with methods of ref- 
erence 2. Methods of applying these results to the calculation of both 
the large and the smaller repeated gust loads in missile operations are 
then developed. 


A 

®n 

b 


C 


SYMBOLS 

gust-response factor, °yj°v 
airplane vertical acceleration, g units 

scale parameter in probability distribution of root-mean- 
square ■ gust velocity 

K g p 0 mSV e 

gust-response factor. 
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c 

D(y) 


Dp 

El, E 
E(<%) 


f (<v) 
0( ) 

Go 

g 

h(q) 


mean geometric chord, ft 

average flight miles required to exceed given values of 
response quantity y 

flight distance in gust-critical flight segment, miles 

functions of flight distance and P ex (defined in equations 
that follow equations (29) and (30), respectively) 

cumulative probability distribution of root-mean- square gust 
velocity 

probability density distribution of root-mean- square gust 
velocity 

average number of peaks of specified response per mile of 
flight exceeding given values of argument 

* •<% 

average number of peaks of specified response per mile of 
flight 

acceleration due to gravity . 

frequency-response function 

gust-response factor 


K 

% 

k 

L 

m 

N 0 

P 


airplane mass ratio, 

girpSc 

gust-response factor (ref. 5) 

turbulence intensity factor describing variations with 
altitude 

scale of turbulence, ft 
slope of lift curve per radian 

average number of peaks of specified response per second of 
flight 

proportion of total flight time or distance in turbulence 
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p ex ( ) probability of exceeding specified value of argument 
g uing area, sq ft 


s = c/L 


U de 

derived gust velocity, fps 

V 

true airspeed, fps 

V e 

equivalent airspeed, Y\jpjp Q 

W 

airplane or missile weight, lb 

y 

response quantity 

y L 

specified value of a response quantity 

Pg 

airplane mass parameter, 

p 

air density, slugs/cu ft 

Po 

air density at sea level, slugs/cu ft 

a an 

root -mean- square normal acceleration 

V 

root -mean- square gust velocity 

a y ■ 

root-mean- square response y 

®(ft) 

power- spectral-density function 

ft- 

frequency, radians /foot 

Subscript 

s : 

1 

nonstorm turbulence 

2 

storm turbulence 
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GENERAL APPROACH 


Turbulence Model 

Tbe approach to be followed in the present study is basically that 
utilized in references 1 and 2. In reference 1, a simplified model was 
used to describe the turbulence experienced in normal airplane operations 
In essence , this model assumed that the turbulence experienced in normal 
operations could be broadly considered to be of two general types: one 

consisting of a severe turbulence condition, represented by turbulence 
encountered in thunderstorms, and termed 11 storm" turbulence and the other 
consisting of a considerably less severe condition, perhaps representa- 
tive of conditions in moderately rough clear air, and termed "nonstorm" 
turbulence. The turbulence for these two conditions was described by 
appropriate average frequency distributions which defined the average 
number of gusts per mile exceeding given values of derived gust veloc- 
ity U^ e . On this basis, the turbulence for a given operation or set 
of atmospheric conditions may be viewed as being given by the following 
relation: 

5 ( U de) = PA( U de) + P 2L( U de) A) 


where 



overall frequency distribution of U^ e encountered in a 
given operation or part of an operation and normally given 
in terms of the average number of gusts per mile of flight 
exceeding a given value 



frequency distribution of U^ e 


for nonstorm turbulence 



frequency distribution of U^ e 


for storm turbulence 


Pi, P^ proportion of total flight distance in nonstorm and storm 
turbulence, respectively 

The appropriate values of Pq_ and P 2 and the appropriate dis- 
tributions of and G^Ude) can conceivably vary with atmos- 

. pheric conditions. Some of the parameters which could affect these 
quantities are 
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Altitude 

Latitude 

Surface conditions (land, water, smooth or rugged terrain) 

Seasons of the year 

Route of operation 

Although efforts to evaluate the variations in turbulence frequency 
and intensity have been made for each of these parameters, no large and 
persistent differences have as yet been established for any of these 
parameters except altitude. For this parameter certain trends appear 
well established, as indicated in a subsequent section. For the other 
parameters the lack of any clear patterns has been, in part, a conse- 
quence of the limitations in the available data which are mostly con- 
fined to operation within and close to the United States. (See, for 
example, refs. 5 and 6.) Also, in many cases the records covered a 
variety of operating conditions in regard to locale, latitude, and even 
seasons of the year, and no separation of the data was possible. In 
several investigations, direct comparisons of turbulence experienced 
at different seasons and on different routes were made and indicated that 
some differences were present. However, the differences observed were 
neither large nor consistent and thus appeared of secondary importance. 

As a consequence of the foregoing limitations in the data, the 
current information on turbulence is restricted to variations in the 
overall turbulence pattern with altitude. Estimates of the quanti- 
ties Pq, P 2 , ^l( u de)> 811(1 ^2( U de) and their variation with altitude 

were given in reference 1 for use in transport -type operations. These 
estimates were based on the limited data available at that time. Since 
that time, a large amount of additional data has been collected, partic- 
ularly for flight altitudes above 10,000 feet and up to altitudes of 
53,000 feet. Also, the data on thunderstorms have since been examined 
in greater detail in reference 7* 

For the foregoing reasons, it appeared appropriate first to revise 
the estimates given earlier in reference 1 for airplane operations. In 
addition, for the present purpose of missile application, adjustments 
to these results are required to account for the effects of storm- 
avoidance procedures normally followed in the airplane operations from 
which the gust data were obtained. 


Power Spectral Representation 

The description of the turbulence in terms of distributions of 
derived gust velocities, as given by equation (l), is then converted 
into a form appropriate for use in power spectral response calculations 
in accordance with the general methods outlined in reference 2. This 
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ccnversiOH provides a turbulence description in terms of the probability 
distributions of the root -mean- square gust velocities. The turbulence 
representation obtained for given flight operations in this manner can 
be expressed in a form analogous to that given by equation (l) as 


where 

f(v) 

p l(°w) 

*2(<V) 


^(°v) “ VlK) + p 2 P 2( a w) 


( 2 ) 


probability density distribution of root -mean- square gust 
velocity 

probability density distribution of root -me an- square gust 
velocity for nonstorm turbulence 

probability density distribution of root -me an- square gust 
velocity for storm turbulence 


As m equation (l), and P 2 represent the proportion of total 

flight time . spent in nonstorm turbulence and in storm turbulence, respec- 
tively This conversion is performed on the basis of an assumed power 
spectral shape as in reference 2. 


Gust Response Calculations 

The representation of the turbulence environment in the form of 
equation (2) can then be applied to the problems of gust response calcu- 
lations by utilizing the general methods described in reference 2 As 
indicated therein, for given conditions the expected response history 
111 7 01 the ^ airplane (where y may be taken as the airplane accel- 

eration, bending moment, stress, or any response quantity) is given by 


G(y) = G f 


r ro - 

Jo 


-y 7 2 v 2 a 2 


ddw- 


( 3 ) 


where 

G(y) 

G^ 


average number of response peaks per mile of flight exceedi 
given values of y 


average 

air 


number of response peaks per mile of flight in roug 
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A ratio of root-mean- square values of specified response y 

and vertical gust velocity (for a given airplane and set 
of conditions and within the framework of linear theory, 
this value depends only upon the form of the gust spectrum), 
°y/ a w 

In the present study, the choice of appropriate functional forms 
for the gust distribution yields a simple form for the proba- 

bility density distribution of root-mean- square, gust velocity f(oy), 
which, in turn, permits a closed-form integration of equation (3) that 
yields a number of results that permit rapid estimation of the large * 
and the repeated gust loads. 


TURBULENCE ENVIRONMENT 


In this section, flight measurements of atmospheric turbulence 
reviewed and a description of the turbulence environment is derived 
terms of the quantities defined in equation (l) ^Pq, P 2 , ^i(% e ), 

G2(Ude)) and in equation (2) (fi(cfy) and f 2 (c%)). 

Nonstorm and Storm Gust Distributions and G 2 ( u de) 

Flight measurements of the gust-velocity distributions indicate that 
the nonstorm and storm gust distributions Gq^Uq^^ and G 2 ^Uqg^, respec- 
tively, vary widely from one day or condition to the next. They do, 
however, on the average show persistent trends with altitude. In refer- 
ence 1, two basic distributions, herein designated by Gq*(u de ) and 

G 2*( U de)> were chosen on the basis of the data available at that time 
and estimates were then made of their variation with altitude. In these 
terms, the distributions G l( U de) and G 2 ( U de) for a given altitude 
are given by 

G i(Ude).= G *(%e/ki) U = l,2) ( 4 ) 

where the quantity k is an intensity parameter which varies with alti- 
tude. The basic distributions Gq^^Ude^ and G 2 *(^Ud6^ used in refer- 
ence' 1 are given in figure 1. The variations in k for the two types 
of turbulence are designated by kq and k 2 and the results used in 
reference 1 for the variations with altitude of these two quantities 


are 
in , 
and 
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are shown in figure 2. Note that for the stom- turbulence case, the 
intensity was taken as the same at all altitudes ( kg = 1.0). 

As a part of the present study, a review was made of the more recent 
data on the variations of turbulence with altitude. This review indi- 
cated that a minor modification in the choice of G-j_* (u^ e ) was desirable 

in order to reflect more closely the values of G^O) (subsequently des- 
ignated as G 0 ) measured in flight tests. The modified distribution 
^l*( U de) is s ^own in figure 1 and is given by 

• G l*( u de)= 20e ' (5) 

The estimates of kq given in reference 1 were, however, retained 
unchanged except that estimates for the lower altitudes (0 to 5,000 feet) 
were added and are shown in figure 2. This extension was made in order 
to represent more adequately conditions at very low altitudes which 
appear of particular interest in certain applications. 

In regard to the distributions of storm turbulence Gg(U de ^ , it 

appeared appropriate to modify the distributions utilized in reference 1, 
as indicated in figure 1, in order to reflect more closely the results 
obtained in reference 7 * The curve shown is based on the results given 
in table III of reference 7 an-d represents a more severe turbulence con- 
dition than that given in reference 1. In addition, this modification 
has the additional advantage for present purposes of yielding a simple 
exponential form for the distributions of Gg*(U de ) (as can be seen from 

the straight-line character of the curve on semilogarithmic paper). The 
distribution is given by 

G 2*( u de) = 15 e'W 5-3 (6) 

The more severe turbulence condition represented by the present choice 
is, by itself, not significant inasmuch as the storm turbulence that 
applies to operations depends also on the values for Pg. 

In addition, the results of reference 7 suggest that for altitudes 
above 20,000 feet the intensity of the turbulence decreases with increase 
in altitude. This result is in accord with the general impression of 
many pilots and is consistent with what may be expected from meteoro- 
logical considerations. (The relatively low moisture content and greater 
stability of the atmosphere at the higher altitudes would tend to make 
smaller amounts of energy available for vertical and turbulent motion.) 
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As a consequence, it appeared reasonable to allow k 2 

altitude above 25,000 feet, as indicated in figure 2. 
however, is arbitrary. 


to decrease with 
The choice made. 


, „ C( ?f^ inin S the results of figures 1 and 2 in accordance with 
equa ion ), the distributions Gi(U^ e ) and G 2 (U^ e ) appropriate for 

each of the altitude brackets are obtained and are shown in figure 3 
1 ° west lc 6 000 feet, separate distributions are shown for the 
altitude brackets of 0 to 2,000 feet and 2,000 to 10,000 feet. These 
frequency distributions are given as follows: 

G l(%e) - 20 e'^/ 2 - 2 ^ (?) 

and 


G 2( u de) ..= .Ue 


-U de /5-3k, 



where the values of k x and k 2 for the various altitudes are defined 
in figure 2. 


, ri I ? 1S interest to note that the coefficients (to be designated 
by G 0 J on the right-hand sides of equations (5) to (8) - namely, 

15 and 20^- define the average number of gust peaks per mile. Except 
or the difference in units, this quantity is approximately related to 
the characteristic frequency N Q of reference 2 (the number of positive 

acceleration peaks per second). These definitions imply that 


° (2) (1.467) (36OO) G ° 10560 G ° ^ 

where V. is the airplane speed in feet per second, and the coefficient 
arises from the fact that N Q is based on positive peaks only, whereas 
Gq and the gust data include both positive and negative peaks. A char- 
acteristic value for 'the airspeed V for the airplanes used in the gust- 
data collections is about 350 feet per second which yields values of N Q 

of about 0.5 and 0.7. These values are reasonably consistent with the 
estimates of N 0 given in reference 2 for most of the airplanes con- 
sidered therein, and, thus, the relation of equation (9) is assumed to 
apply m subsequent considerations. 
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The foregoing estimates of the gust distribution were, in large 
part,, based on center- of -gravity normal-acceleration measurements obtained 
from transport operations. It is well known that for many airplanes 
the .effects of airplane flexibility give rise to substantial amplifi- 
cations of the airplane center- of -gravity accelerations. As a conse- 
quence, the values of U de derived from such amplified accelerations 

would likewise tend to be amplified. In reference 1, a simple correction 
or reduction of 20 percent was applied to the acceleration measurements 
and thus to the gust velocities to account for this effect. In the 
present investigation^ the same correction was used in the determination 
of the distribution Gq(u de ) • However, for the distribution 0 2 (u de ), 

no such correction was believed necessary inasmuch as the airplanes used 
in obtaining most of the thunderstorm gust data were relatively stiff 
and dynamic effects on the center-of-gravity accelerations were small. 

In comparing these distributions with operational data, this difference 
must be kept in mind and the effects of flexibility on the operational 
data be considered. " — 


Proportions of Flight Distance in Nonstorm and Storm 
Turbulence P d and P 2 

In order to determine appropriate proportions of flight distance in 
nonstorm and storm turbulence P d and P 2 for transport operations, 

equation (l) was used with the results of figure 3 to approximate the 
gust distributions measured in transport operations. Simple graphical 
procedures were used and yielded estimates of P d and P 2 which gave 

good, representations of the measured data. Inasmuch as the data from 
various operations for a given altitude bracket varied widely, average 
values of P d and P 2 were obtained. The' values of P-]_ and P 2 

obtained for the various altitudes are shown in figure k. For compari- 
son, the values of P d and P 2 from reference 1 are also shown. The 
same 20-percent correction, discussed previously, to account for dynamic 
WQ’S also applied to the operational gust data in deriving esti- 
mates of P d and P 2 . 

Inasmuch as the operational data available for the higher altitudes 
(above 20,000 feet) were limited, estimates of P d and particularly 

of P 2 are at best crude. In estimating values of P 2 , no flight data 

were available and recourse to indirect evidence such as that given in 
reference 8 on the distribution of thunderstorm cloud tops was necessary. 
These data were used as a basis for extrapolating the values of P 2 

obtained from the gust data for the lower altitudes to the higher 
altitudes. 
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The present values of P d and Pg shown in figure 4 differ in a 

number of respects from those given in reference 1. In regard to the 
values of P d , the most significant difference is the increase in value 

for the altitudes between 20,000 and 40,000 feet. This increase is 
indicated "by recent unpublished studies and in reference 9 and is asso- 
ciated with the presence of the jet stream in this altitude range. The 
decrease in values of Pg from those given in reference 1 also appears 

large. However, the distribution Gg(U de ) has, for present purposes, 

been selected to be more severe than that used in reference 1. The net 
effect of these two modifications is a small increase in the severe tur- 
bulence condition for the present estimates. 


Overall Gust Distribution G(U de ) 


By combining the results obtained in figures 3 and 4 in accordance 
with equation (l), the overall distributions of gust velocity G(U de ) 

for the various altitude brackets are obtained and are given in figure 5. 
For this purpose, average values of P d and Pg for the various alti- ■* 

tude brackets were determined from figure 4 . The actual values used 
are summarized in the following table: 


Altitude, ft 

p i 

P 2 

0 to 2,000 . . . 

0.32 

0.00025 

2,000 to 10,000 . . . 

0.08 

0.0008 

10,000 to 20,000 . . . 

0.045 

0.0004 

20,000 to 30,000 . . . 

0.06 

0.00013 

30,000 to 40,000 . . . 

0.065 

0.000045 

40,000 to 50,000 . . . 

0.023 

0 . 00001 

50,000 to 60,000 . . . 

0.02 

0 


The frequency distributions of figure 5 are all given by the following 
expression: 


G( U de) 


20Ppe 


de 


A 


-U, /2.2k 


k l 


+ 15P 2 e 


; u d e/5- 


3K 


( 10 ) 


where the values of P d and Pg are those given in the foregoing table 
and the values of k d and kg are obtained from curves in figure 2 at 
the 'midpoints of the various altitude brackets. 
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Modifications To Account for Storm-Avoidance Procedures 


For purposes of missile applications, modifications are required to 
the foregoing results in order to elimiriate the effects on the data of 
airplane storm-avoidance procedures. These modifications can at best be 
only crudely estimated on the basis of available information. Available 
information indicates that little effective effort is nor mally made by 
pilots to avoid the lighter or nonstorm-turbulence areas. However, serious 

Sttle r au^tit C at 1Ve eff 7 tS ?" £ n0mally mde t0 aVOld atom-turbul^ce areas 
Little quantitative information is available on the consequences on the gust 

thnt r h!e Ce ° f auch . stor m-avoidance procedures. Some indirect information 
that has some bearing on this problem is, however, available and includes 
data on the frequency of thunderstorms, their average horizontal dimensions 

tha^thun d^ a + 10nS ’ thSir altltude extent. Roughly it is estimated 

United States a^Tb° CCUr ’ °“ aVerage ’ °n about 30 days per year for the 
United States and have an average duration of perhaps two hours. It would 

thus appear that for a given location the probability of a thunderstorm 


being- present is approximately equal to - (S 0 )( 2 ) 


or O.OO7. Comparison 


(360) (24) 

of this value with those of figure 4 for airplane operations suggests that 
airplanes may veil avoid a large part of the atmospheric storms. Inasmuch 
as thunderstorms are probably less frequent on a vorldvide basis, somevhat 
lover values than 0.007 appeared appropriate for present purposes. The 
values of P 2 selected as representative for missiles in all-veather 
operations are those shown in figure 4 (b) . 

Application of these modified values of P £ in equation (10) yields 
the distribution G(u de ) appropriate for all-weather missile operations, 

and these distributions are given in figure 6. In general, they repre- 
sent a more severe gust history than that given earlier for airplane 

operations and for the less frequent gusts, say, = 10"5, are 

roughly 40 percent more severe at the various altitude levels . Analyt- 

eive?I»ri Se d ! strlbu y° ns “ay be represented by the same expression as 
given earlier m equation (10). 


Conversion to Pover Spectral Form 

The distribution cfu de j in figure 6 may be converted into a form 

appropriate for power spectral response calculations by making use of 

^oli a o? r ?he t° f b r f erence 2 ’ As lndicated therein, if the power spectral 
° d dhe turbulence is assumed invariant, the turbulence history expe- 

bition Oft f 1 alrd,lane “y be defined by the probability density distri- 
bution of the root-mean-square gust velocity f(p). On the basis of 
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the results of reference 2 and as given in equations (3)_and ( 9 ) herein, 
f is related to the airplane acceleration history G(a n ), in terms 

of the average number of acceleration peaks per mile exceeding given 
values of a^, by the relation 


G(a n ) = 


10560 


f (cv)e 


A 2 


N Q characteristic frequency of airplane acceleration response 

and approximately specifies average number of peak accel- 
erations per second 

A = Oan/°w for the specified airplane and spectral form 
Inasmuch as 


t _ 2W n _ 1 

de K-p mSV e n C n 
go *= 


_ K gPo mSV e 

where C - — — ■, the derived gust velocity may be viewed as a 

reduced or normalized acceleration and the distribution of peak values 
of U^ e is, in turn, from equations (ll) and (12) given by 


G ( U de) = 


„ -Ude 2 /2o w 2 (A) 2 

? W e I 

v j 0 


From equations (lO) and ( 13 ), the distributions f(a w ) and. G(% e ) are 
related by 


-U,q 2 


' \C J 


! 056 on 0 


-U de 2.2k x 


20P x e 


!5P 2 € 


-%e/5-3k 2 
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where , as indicated earlier (eg. ( 9 )), 


V 

10560N o 


— and — 
15 20 


for the airplanes used in the gust measurements. Thus, to this 
approximation 



The solution of equation ( 15 ) is given hy 





f (a v ) - Pi?i(%) + P 2 f 2 ( CP w) 

where 


(15) 


(16) 





1 6 


MCA TN ^332 


The determination of the values of 


thus depends upon 


which is the ratio of the acceleration response to unit discrete gusts 
of the standard form (cosine shape) and the root-mean-square acceleration 
response to a random gust input a v = 1. This ratio must be established 

for the airplane involved in the gust-data collection program. 

For the single-degree-of -freedom case, vertical motion only (which 
appears adequate for present purposes as indicated in ref. 2), 


)VSm fifes) 


where the term 


fifes") 


is a gust-response factor depending on K, the 


airplane mass ratio, and on s, the ratio of wing chord to scale of tur- 
bulence L. (See refs. 2 and 10.) Thus, from equations (l2) and (17), 


fifes)/* 


For present purposes, this ratio was evaluated on the basis of a char- 
acteristic transport configuration as given in table V of reference 2 
in order to determine values of b^ and of the Northrop P-61C airplane 

(the airplane actually used in the Thunderstorm Project gust survey) for 
the determination of bg. The same form of gust power spectrum as that 

in reference 2 was used as well as a value of the scale of turbulence L 

Q 

of 1,000 feet. The ratio — varies with altitude and the actual values 

A 

obtained are given in table I. The values of b 1 and b 2 for the var- 
ious altitude brackets are also given in the table. The associated 
probability density and cumulative probability distributions f (<v) and 

F fe) for the various altitude brackets are given in figure 7. The dis- 

• • ^ / % A 

tributions of a w for the nonstorm turbulence fq fe) and FqfcO and 
the storm turbulence fg^w) -^2 (°w) 3X6 a I so given separately in 

figure 8 for each of the altitude brackets. 
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application to gust- load calculations 


In the preceding section, a simplified description of the turbulence 
at the various altitudes was derived in terms of the probability density 
distributions of the root-mean- square gust velocity. This distribution 
is given by 


f (^w) - P 1 




+ P 2 



( 19 ) 


in which the parameters P^ and P 2 represent the proportion of flight 

time (or distance) in nonstorm and storm turbulence, respectively, and 
bp and b 2 represent scale-parameter values for the individual prob- 
ability distributions of o v for the two types of turbulence. The 
values of P^, P 2 , b^ and b 2 varied with altitude. In this section, 

the foregoing specification of the turbulence environment is applied to 
the problems of missile gust- load-history calculations. 


Estimation of Severe Gust Loads 

As indicated in equation ( 3 ), the gust response history for a given 
airplane under given conditions, exposed to a gust history consisting of 
a series of locally stationary Gaussian processes of common spectral form 
(as, for example, defined by eq. ( 19 ) ) ^ may in general be expressed as 


G(y) - G, 


, GO ^ __*\r2 

f (<V)e 


-y^/2i% 2 A 2 




( 20 ) 


where 

y 


response quantity of concern (load, bending moment, stress, 
and so forth) 

number of response peaks per mile of flight in rough air 


A = a. 


r y/ a w for ' t ^ ie specified spectral form of the gust input and, as 

indicated in reference 2, need not be restricted to single- 
degree-of -freedom systems 
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Substituting equation (19) into equation (20) and integrating yields 

-y/biA -y/tpA 

<J(y) = p l G 0 e + p 2 G o e (21) 

The results given by equation (21) may be viewed as a description of 
the statistics of the peak values of y and represent averages for 
extended operations under the specified conditions. As such, they do 
not apply directly to a single missile flight "but must "be viewed as the 
overall response histories of a large number of missiles for the specified 
conditions. 

Equation (21) must be applied separately to each significant segment 
of the flight plan since the turbulence parameters P^> P 2 , b 1 , and b2 

vary with altitude and the missile parameters G Q and A may also be 

expected to . vary widely with the flight segment. If several flight seg- 
ments are significant, either the overall load history G(y) must be 
determined as a weighted average (weighted, perhaps best, by the flight 
distances in each segment) or the load histories for individual flight 
segments must be considered separately. In many practical cases, one 
or two flight segments only are gust critical. This condition simplifies 
matters appreciably and is considered in a subsequent section. 

If the load history, as specified by equation (21), is examined, 
it is clear that a small but finite probability of exceeding large values 
of y exists no matter what values of y are chosen. In any case, it 
is therefore impossible to select a value which will never be exceeded. 
Instead, it is necessary to accept some tolerable risk level or some 
finite probability of exceeding a chosen value. The actual probability 
value chosen would presumably depend upon the particular missile, the 
consequences of a structural failure, and economic and military tactical 
considerations . The question of the choice of the probability value Is 
beyond the scope of this paper, and consideration herein is restricted 
to the problem of determining the load value once the probability of 
exceedance is chosen. 

Consider the case of a single missile flight involving a flight 
distance Dp. This flight may be viewed as yielding a sample of the 
random process y(t) of distance Dp. The random process y(t) has 
an average of one exceedance of a specified value y^ in D (yj) flight 
miles where 

5 M • <-> 
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and where G(y) specifies the average load history of the missiles for 
extended flights. If it is assumed that the exceedances of y^ are 

distributed at random, then the probability of exceeding y^ in a given 
flight distance D r is approximately given by 


(23) 

provided that b(y L ) » D p which is assumed to be the case of interest. 
(The assumption of random distributions of y^ on y(t) does not apply 

in a strict sense to the random process y(t) as specified by the non- 
stationary input of equation (19)* This assumption and the approximation 
of equation (23) are adequate for present purposes and are conservative 
to the extent that the cases of multiple values of y T separated by a 

flight distance less than M are excluded. ) 

For given values of D r and P ex (y L ), equations (22) and (23) 
specify the value of The result of the load calculation given 

"by equation (21) may then be used to determine the required value of y 
to achieve the desired ^ex^L)* several flight segments are being 
evaluated separately, the value of y^ may be determined in such a man- 
ner that the desired exceedance rate ^exC^l) is given by 



(2k) 


MM, “ “ 1 " ai ’' ia ” 1 

flight segments 'and the probabilities in the various segments are 
assumed independent. 


A Simple Formula for Estimating Severe Gust Loads 

In many cases of interest only a portion of the flight path or a 
single flight segment may be gust critical. If only a single flight 
'segment is gust critical, it appears possible to derive a relatively 
simple formula for y^ in terms of a few significant quantities. For 
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this purpose, it is of interest to examine the relative contributions 
to G(y) of the two terns on the right-hand, side of equation (21). Let 


G(y) 


= Gi(y) + CL>(y) 


(25) 


where 


G x (y) = Pie 


-y/ b i A 


r 


G 2 (y) = ?2 e 


-y/v 


The gust data presented earlier indicate that for the significant alti- 
tude brackets 


P 1 ** 20P 2 


b 2 - 3b! 


(26) 


For these conditions, the relative contributions of the two terms are 
schematically illustrated by the following sketch (a logarithmic scale 
applies to the ordinate): 
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As can be seen from the sketch, the principal contribution to G(y) 
arises from Gq(y) (the nonstorm-turbulence contribution) at low values 
of X and from G 2 (y) (the storm-turbulence contribution) at high values 

xi 

• Which of these two cases is of concern would appear to depend, 

in large part, on the particular missile and the desired exceedance 
rate. It is believed that the region of high values of X is of prin- 

A 

cipal concern although, in some applications where operational consid- 
erations permit planning for the avoidance of storm turbulence, the 
Gp(y) case may alone be applicable. 

In either case, equation (21) yields 


y L = t>±A log e + log e G 0 - log e G(y L )j 
yL = jlogg P i + logg G c + log e D 


(i = 1,2) (27) 


Substituting for D^y L j from equation (23) into equations (27) yields 


XL = b i A lo G e P i + lo S e G o + lo S e “ (28) 

P ex(^p) 

which is a simple and useful result. Equation (28) specifies a value 
of y^ in terms of the following groups of parameters : 

(a) Gust input parameters bj_ and 

(b) Missile response dynamics A and G 0 

(c) Operational parameter D r 

(d) Desired exceedance rate P QV 

CA 

From figure 4 and table I, representative values of P and b for 
the altitude brackets of 0 to 40, 000 feet are for the nonstorm-turbulence 
case 

P x = 0.06 b-j^ = 3.5 
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and for the storm- turbulence case 


P 2 = 0.0025 


b 2 = 10.5 


■I 


Utilizing these values in equation (28) yields the following results: 
For the nonstorm- turbulence case, 


where 


y L = 3.5A|log e 0.06 + log e G q + log e 

\ ^ex 


= 3-5A^log 


e G o + E l) 


( 29 ) 


D, 


E 1 = l°g e - 2.81 


ex 


m 


For the storm- turbulence case. 


D, 


y = 10. 5 A log e 0.0025 + log e 0 o + log e 


v 


y L = 10 . 5A (log e G 0 + 


ex 


( 50 ) 


where 


E 2 = log e 


D 


ex 


6.0 


The values of E^ and E 2 are shown in figure 9 for a range of values 
of> P ex (^l) from 0.001 to 0.2 and for a range of values of Dp from 10 

to 5; 000 miles. The charts of figure 9 can be used directly along with 
the missile response parameters A and G 0 to determine the load values 
in accordance with equations (29) and (30)- The simple form of these 
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results suggests that they could he used in preliminary design studies 
and in the development of design specifications. 

In order to illustrate the applications of the foregoing results 
an example is given. Consider a missile having a flight range in the' 
lower atmosphere D r of 100 miles and values for G 0 of 10 and for 

P ex (^ l) of °- 01 - For this case, a value of E 2 of 3.2 is obtained 

from figure 9(h). Using these values in equation (30) for the storm- 
turbulence case yields 


y L = 10.5A(log e G q + 3.2 ) = 58a ( 31) 

It is of interest to note that doubling the range D r to 200 miles 
yields 


y L = 64 a 

or about a 10-percent increase in the value of y L# (a 10-percent 
increase is also obtained if P ex (y L ) is reduced by one-half, that is, 

Pexta) - 0-005.) 

If the missile operations are restricted to the avoidance of storm 
areas and equation (29) for the nonstorm-turbulence case is considered 
applicable, the value obtained for the initial example is as follows: 


y L - 30A 


It^is clear that a large reduction in the value of y^ (from 58A to 

30A) may be achieved by the avoidance of storm- turbulence areas. The 
structural penalty for all-weather missile operations thus appears large. 


Estimation of Repeated Gust Loads 


• P ro ^ lem of* calculating the repeated loads and developing a 

atigue loading differs in a significant respect from that of the limit 
l Q ad case. In the case of large loads, it is useful to consider the 
verall history of a fleet of missiles to insure that, on the average 
e critica l load is exceeded with a given frequency. In the fatigue' 
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case, the fleet concept in this form cannot he used. Instead, the 
cumulative load history of the individual missiles is of concern. The 
determination of such cumulative load histories requires information 
on the concurrent gust histories for the various flight segments of a 
particular missile flight. No information of this type is available. 

In some practical cases, significant simplifications may he feasible. 

One such possible simplification is considered herein. 

It is assumed that the missile gust history for the significant 
part of the flight is statistically homogeneous and is specified by 
a given value of the root-mean-square gust velocity. This assumption 
may be expected to apply best to the case of missiles of short flight 
duration and appears, in general, to be conservative. On this basis, 
the cumulative load history for a given missile may be obtained from the 
following equation: 


Q-t(y) = ) DiGi(y) 


(32) 


where 

Gfc(y-) 


expected number of response peaks exceeding given values of y 


Dp flight distance in ith flight segment 


Gi(y) 


response history in ith segment which is obtained from 



This procedure assumes that the root-mean- square gust velocity is con- 
stant but that G Q and A vary with flight segment, (it also assumes 

that the flight distance is suff iciently_long to insure that the load 
history is close to the expected value Gt(y) • ) For a given missile 
operation, the loqd history (and thus the fatigue damage) from equa- 
tion (32) depends only upon a w . The distribution of the load histories 

for a series of missiles, in turn, depends upon the probability distri- 
bution of o^. Thus, the specification of a value of a v which is 

exceeded with a given desired probability implies that the associated 
load history, as given by equation ( 32 ), is likewise exceeded with this 
same probability. For example, for a probability level of 0.001, fig- 
ure 7 (b) indicates that the value of a w exceeded with this probability 

varies between 6 and 11 for the various altitude brackets (ignoring the 
lowest altitude level). The conservative choice of a value for c v of 
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11 feet per second for calculations of repeated loads in equation (32) 
would thus yield a load history which would he exceeded with a proba- 


bility of less than 0.001. 


COMMENTS ON APPLICATIONS AND LIMITATIONS 


The applications of the results obtained in the previous section 
to load calculations pose a number of problems. The applications, in 
general, require the determination of the quantities A, G Q , D r , 

and P ex . The choice of values for the last quantity P ex depends upon 

the particular problem and need not be of concern herein. The remaining 
quantities A and G 0 , which define the missile response characteristics, 
and D r , which depends upon the operational flight path, warrant some 
comment . 

The quantities A and G Q , in practice, probably have to be deter- 
mined by analytic means although, in some cases, direct experimental , 
determinations may be possible. Analytically, these quantities may be 
defined as follows (ref. 2): 



r poo 0 ql/2 

G ° = iy J 0 dn (3^+ ) 

where 

$ w (2) power spectrum of gust velocity 

H(fi) frequency-response function of missile, defining specified 

response of missile to unit sinusoidal gusts of frequency 2 

As specified by equations ( 33 ) and ( 34 ), no limitations exist, other than 
the usual one of a linear system, in the determination of H(G). Thus, 
in addition to the translational and rotational degrees of freedom, the 
effects of the missile control system and structural dynamics may be 
included in the analysis. 
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The determination of the appropriate value of D r for a given 

missile operation may, in practice, also involve some difficulty. As 
utilized in equations (29) and (30), D r is the flight distance in the 

gust-critical flight section. Inasmuch as the turbulence decreases 
rapidly above 40,000 to 50*000 feet, an upper limit in the value of -D^ 

is the total flight distance below, say, 50,000 feet. In addition, only 
a small part of this flight distance may be at relatively high dynamic 
pressure. As a consequence, some arbitrary criterion for the determina- 
tion of 1^, such as the flight distance below 50,000 feet and within 
20 percent of the maximum dynamic pressure, may be desirable. 

Although principal consideration has been given herein to the prob- 
lems of gust-load calculations for missiles, the present results may 
also find application to other problems, such as the esti ma tion of 
missile-motion response histories which may be required in guidance and 
tracking studies. In addition, many of the present results can, with 
minor modifications, be applied to airplane operations. For example, 
for transport-airplane operations (without radar for storm- turbulence 
avoidance) the appropriate value of P 2 in equation (30) differs from ' 

the value used for the all-weather missile case and instead would be 
based on the values given for airplane operation in figure 4. 

The foregoing analysis based on turbulence data collected by air- 
planes in horizontal flight applies best to the case of missiles in 
flight paths similar to those of airplanes - that is, flight operations 
involving horizontal or moderately inclined flight paths. However, a 
large number of missile missions require rapid exit and entry through 
the lower atmospheric layers where air-motion disturbances are likely 
to give rise to significant loads. Missiles in .such flight operations 
are likely to have near- vertical flight paths. For these cases, the 
use of airplane data is open to question for several reasons. First, 
the assumption of even local isotropy is probably most closely approxi- 
mated in the atmosphere for horizontal layers and is unlikely to apply 
very well to the case of vertical flight paths because of the rapid 
changes in mean wind flow with altitude. This is particularly evident 
when it is recalled that turbulent areas are normally layers with a 
horizontal extent of 10 to 100 miles and with relatively thin vertical 
thicknesses of only several thousand feet. In addition, the rapid vari- 
ations in horizontal wind speed with height (sometimes reaching values 
of 100 miles per hour in a few thousand feet as in jet-stream areas) 
are of an order of magnitude larger than the vertical gust velocities 
encountered in horizontal flight. These large wind shears exist at 
altitudes of 5 to 10 miles and appear to be the principal source of 
atmospheric disturbances applicable to missiles in vertical flight. 

In addition to those difficulties, missiles in vertical flight normally 
undergo such rapid variations in airspeed, dynamic pressure, and air 
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density that it is questionable whether a locally time -invariant -system 
approach, as utilized herein, would apply. For these reasons, it is 
felt that the case of missile operations in near -vertical flight paths 
requires a separate and different approach centered upon direct meas- 
urements of the variations in horizontal wind with altitude as distinct 
from the measurements of turbulence obtained from airplanes in horizontal 
flight. 

CONCLUDING REMARKS 


Airplane measurements of atmospheric turbulence have been utilized 
to derive a simplified description of the atmospheric turbulence environ- 
ment appropriate for missile operations. This description was then 
applied in developing an approach to the estimation of severe and repeated 
gust loads. Relations are given for calculating severe loads that are 
exceeded with a given probability as a function of turbulence parameters, 
the missile response characteristics, and the flight distance. Results 
are given for two cases: one which might be considered an all-weather 

operation and the other a limited-weather operation involving the avoid- 
ance of storm-turbulence areas. The levels of load values obtained for 
the two cases differ by a large amount. A simple procedure for esti- 
mating the repeated gust-load histories for missiles is also given. 

Inasmuch as the present results are based on airplane measurements 
obtained in essentially horizontal flight, they appear applicable to 
missile flight operations involving only horizontal or near-horizontal 
flight. They do not, in particular, appear well suited for missile 
operations involving near-vertical flight paths through the lower atmos- 
phere. For such operations, the changes in the horizontal wind with 
altitude appear to be the largest source of air -mot ion disturbance. 

This case appears to require a separate and different approach and one 
based on direct and detailed wind-shear measurements. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 1 6 , 1958. 




28 


NACA TN 4332 


REFERENCES 


1. McDougal, Robert L., Coleman, Thomas L . , and Smithy Philip L. : The 

Variation of Atmospheric Turbulence With Altitude and Its Effect 
on Airplane Gust Loads. MCA RM 'L53G15a, 1953. 

2. Press, Harry, Meadows, May T., and Hadlock, Ivan: A Reevaluation of 

Data on Atmospheric Turbulence and Airplane Gust Loads for Appli- 
cation in Spectral Calculations. MCA Rep. 1272, 1956 . (Supersedes 
MCA TN 33^2 by Press, Meadows, and Hadlock and TN 3540 by Press 
and Meadows . ) 

3. Notess, Charles B.: A Study of the Nature of Atmospheric Turbulence 

Based Upon Flight Measurement of the Gust Velocity Components. 

WADC Tech. Rep. 57-259, ASTIA Doc. No. AD 118242, U. S. Air Force, 
May 1957. 

4. Press, Harry: Atmospheric Turbulence Environment With Special Refer- 

ence to Continuous Turbulence. Rep. 115 , AGARD, North Atlantic 
Treaty Organization (Paris), Apr. -May 1957. 

5 . Pratt, Kermit G., and Walker, Walter G. : A Revised Gust -Load Formula 

and a Re-Evaluation of V-G Data Taken on Civil Transport Airplanes 
From 1933 to 1950. MCA Rep. 1206, 1954. (Supersedes NACA TN's 
2964 by Kermit G. Pratt and 3041 by Walter G. Walker.) 

1 

6 . Copp, Martin R., and Coleman, Thomas L.: An Analysis of Acceleration, 

Airspeed, and Gust-Velocity Data From One Type of Four -Engine Trans- 
port Airplane Operated Over Two Domestic Routes. MCA TN 3475, 

1955. 

7 . Tolefson, H. B.: Summary of Derived Gust Velocities Obtained From 

Measurements Within Thunderstorms. MCA Rep. 1285, 1956. (Super- 
sedes NACA TN 3538.) 

8 . Anon.: The Thunderstorm. U. S. Dept, of Commerce, June 1949. 

9. Coleman, Thomas L., and Coe, Emilie C.: Airplane Measurements of 

Atmospheric Turbulence for Altitudes Between 20,000 and 55,000 Feet 
Over the Western Part of the United States. NACA RM L57G02, 1957. 

Fung, Y. C.: Statistical Aspects of Dynamic Loads. Jour. Aero. Sci., 

vol. 20, no. 5 , May 1953, pp. 317-330. 








10 . 



5661 m > vo m 


I 


6*51 

51* 6 

668*5 

80*1 * 

561* 

968* 

0*009 

1’TT 

69*5 

9l5*5 

9I*T 

551* 

018* 

6*516 

T*TT. 

56*5 

161*1 

65*1 

969* 

158* 

1*165 

5*11 

51*5 

161*1 

5^*1 

196* 

861* 

1*091 

1*01 

16*1 

195*1 

96*1 

561* 

89^* 

0*111 

8*6 

68*1 

160*1 

69*1 

661* 

96I* 

6*68 

1*6 

8l*T 

610*1 

91*1 

601*0 

811*0 

6*5l 


1*66i 


000^09 oq 000 'o6 
000 '06 oq 000*01 
000*01 oq 000*06 
000 *06 01 000*05 
000 * 05 oq 000 * 01 
000 *01 oq- 000 '5 
000*5 oq. 0 


* K a 


t d/ ° d r (e'xyi f 9 * WW\ §x x ' §Ti 

[qj 000*1 = 1 lu’BfP'Bj:- .red’ £g*i = m iqj 6*01 = 2 *qj 1*599 = S *01 000*06 = A] 

5q j:aq.aurBXBd sqBog (q) 


8*5 

80*6 

668*5 

10*1 

008* 

668* 

0*666 

o*oli 

6* 

1*6 

66*5 

8l5*5 

11*1 

861* 

018* 

0*566 

6*601 

5*1 

1*6 

15*5 

161*1 

65*1 

099* 

658* 

1*905 

6*69 

6*l" 

6*6 

60*5 

161*1 

96*1 

986* 

861* 

0*611 

1*61 

1*1 

1*6 

98 *T 

195*1 

11*1 

616* 

861* 

8*101 

6*56 

0*5 

8*6 

9l*1 

160 * I 

09*1 

061 * 

55l* 

6*9! 

6*15 

5*5 

9 ,+ f 

11*1 

610*1 

69*1 

651*0 

601*0 

! 

8*69 

0*15 

9*5 

T d 

V 

"S’ 

Q. 

O 

<a^ 


Trip!' 

s x 

X 

»Tl 

^5*5 


000*09 oq ooo *06 

000 *06 oq 000*01 
000*01 oq 000*06 
000 *06 oq 000*05 
000*05 oq 000*01 
000*01 oq 000*5 
000*5 oq 0 


I 1 1 1 1 L— » 1 1 zr 

[qj 000*1 = *1 itrBiP'e.i j;sd 66*1 = m iqj 1*£ 1 = 2 691*1 - S *01 000*11 = ttj 

Iq .isqaurBjred Gxeog (b) 


SMOiinaiiusia iiinavao aa isno mi 


saaiawvava aivos hhi mo MOiivmvAa - * i -aiffVL 






ion: 






'y[ pire sjcaq-surejced souainqjriq. jo spnp'iq.p'e U0Tq.' B J* rB A ~'Z 3-xn^T^ 


q.39j 000'9 9 13 Atsus^ui q.sno/iCq.Tsue^u-p q.sn .9 

0*1 9’ 9* T 




Average- number of gusts per mile, G^U^e 




MCA TN 4332 



Derived gust velocity, U de , fps 

Figure 3. -^Distributions of derived gust velocity for nonstorm turbu- 
lence Gi^U^e) and- storm turbulence at various altitudes. 
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Figure 5.- Distribution of overall derived gust velocity 

airplanes at. various altitudes. S 
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(a) Probability density for nonstorm turbulence f-^ 

Figure 8.- Probability density and cumulative probability d 
of root -mean- square gust velocity for nonstorm turbulenc 
and Fj^cr^ and for storm turbulence 811(1 F 2 (cr 

various altitudes. 
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Figure 8.- Continued. 
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(c) Cumulative probability for nonstorm turbulence 

Figure 8.- Continued. 
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